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Abstract 

"sj- ' The 0(a s 2 ) contribution to the Energy-Energy Correlation function (EEC) jj], §, |[ [jj of 

e + e~ — > hadrons is calculated to high precision and the results are shown to be larger than 
previously reported Ji], [n], [n], |l^, [l3| . The consistency with the leading logarithm approximation 
and the accurate cancellation of infrared singularities exhibited by the new calculation suggest 
that it is reliable. We offer evidence that the source of the disagreement with previous results lies 

^—^ ' in the regulation of double singularities. 

in ■ 

The energy-energy correlation function (EEC) p], |2|, ||, ^| for e + e~ annihilation into hadrons is 
' widely used as a measure of the strong coupling constant a s @, J?J and is potentially one of the 

most precise and detailed experimental tests of QCD available M, S. However, that potential has 
not been realized due to disagreement over the predicted value of the next-to-leading order correction 
in the strong coupling constant ||, [l0|, O, [12], OJ. We report on a new calculation of the 0(a s 2 ) 
term using subtraction for control of infrared singularities. Accuracy was checked at every stage by 
symbolic computation, high precision arithmetic, and human calculation. The detailed cancellation of 
singularities in the complicated four-parton states was carefully tested. A more complete description 
will be presented elsewhere . 

The EEC was invented to take advantage of the asymptotic freedom of QCD by viewing the 
products of c + c~ annihilation with a weighting that favored the most energetic hadrons 0, [| |J. 
Conservation of energy requires all energy carried by quarks and gluons to be transferred to detectable 
hadrons, hence the EEC is experimentally and theoretically defined as 
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where a is the total cross section for e + e~ — > hadrons, E n and p n are the energy and momentum 
of particle n, and E to tai is the center of mass energy of the system. The EEC is free of collinear 
singularities since all parallel momenta are linearly summed flq| . 
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Figure 1: The 0(a s 2 ) contribution to the Energy-Energy Correlation function. For comparison we 
display our results (solid circles), the results of Kunszt and Nason E3| (open squares), and the results 
of Richards, Stirling, and Ellis JTo| (open triangles). B values shown are for five active quark flavors 
or Tr = | (see equation ||). 



After factoring out the trivial dependence on the total cross section and sin 2 \ 
the following pcrturbative expansion in the region < x < 7Tj 



L2 , the EEC has 
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Here (Tq is the leading order total cross section, fj, is the renormalization scale, and /3q is the leading 



coefficient of the (3 function: /3q 



li 



tTr. For QCD in this notation, Cp 



Ca = 3, and 



Tr — hNjr, where Np is the number of active quark flavors at energy E to tai- Analytic calculation of 
A yields § 



Ax) = Cf (1 + w) (1 + 3a;) [(2-6cj 2 )log(l + cj- 1 ) +6^-3] 



(3) 



where w = cot 2 (x/2)- No such analytic expression is possible for B(x)- At 0(a s 2 ), the EEC receives 
contributions from four-parton final states at tree level and from three-parton final states with a 
virtual parton forming one internal loop. The three-parton final states pose little challenge, but the 
integrals corresponding to four-parton states with an external angle fixed at x demand numerical as 
well as analytic calculation. 

To calculate contributions near soft or collinear poles, the four-parton expressions were simplified 
to allow analytic integration in the presence of an infrared regulator e (dimension D = 4 — 2e). Using 
the subtraction method of infrared regulation, the simplified expressions were subtracted from exact 
expressions and the finite difference was numerically integrated without infrared regulation (e = 0). 
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Table 1: The coefficients of the leading log expansion of the EEC at large (B^) and small (B~) 
angles. The expansion is as shown in equation |J. Listed are the exact leading log coefficients and the 
coefficients producing the best fit to Clay and Ellis as well as Richards, Stirling, and Ellis JTtJ . 



Analytic integrals of the three-parton and simplified four-parton expressions (at finite e) were then 
added and the sum was shown to remain finite in the limit e — > 0. As in all previous calculations of B, 
we used the expressions derived by Ellis, Ross, and Terrano (ERT) |l(| for the exact three-parton and 
four-parton final states, but we did not use the ERT simplifications or analytic integrals for reasons 
of maximizing numerical convergence. 

Our results (Clay and Ellis or CE) are plotted in Figure [l] along with the results previously 
reported by Richards, Stirling, and Ellis (RSE) [l(J and Kunszt and Nason (KN) fl^l . The mean 
relative numerical uncertainty in our calculation is 0.3%, while for KN it is roughly 4%, both arising 
from the precision of numerical integrations. This uncertainty is insufficient to explain the roughly 
15% overall difference between KN and CE. While it is possible for systematic differences such as 
these to arise from purely numerical errors, we believe there is an analytic error at the heart of the 
disagreement. 

The only known test of the analytic behavior of B is a comparison with the predictions of the leading 
logarithm approximation for large and small angles H. To determine asymptotic behavior, B(x) was 
calculated over the range |cos(x)| < (1 — 10~ 6 ), and the results were compared to an expansion of the 
form 

km B{ X ) = C F Y j Bf Ml/^)} 3 (4) 

3=0 

where 77 * = ^(l±cos(x)). The coefficients B^ that best fit our calculation were found using an 
unconstrained least squares fit and are displayed in Table |l| (we find that B% — B^ = 0, as expected). 
For comparison, we also show the coefficients derived by RSE (l7[] who reported some inconsistency 
with the leading logarithm approximation. No inconsistency is evident in our data. The previously 
unpublished exact values for B$ are based on our conjecture that the form factor for the EEC is the 
same as that for the second energy moment of the Drell- Yan cross section |l|, . The form factor is 
convoluted with a known parton evolution function pof to produce Bq . 

The discrepancy over the value of Bq is significant. With N F = 4, RSE extracted a value of Bq 
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equal to 29.9 ± 2.9, while our calculation predicts a value of 47.8 ± 0.8 (see Table [j]). Based on our 
preliminary analysis of data from KN as well as Glover and Sutton (GS) we conclude that neither 
is consistent with the values of Bq from either CE or RSE. It is unfortunate that the coefficient that 
best discriminates between the various calculations is unknown. An independent calculation of Bq 
would be very useful for resolving the disagreement. 

To explore the source of the disagreement, we parameterize B as a sum of three functions 

B( X ) = C F ( C F BcAx) + Ca BcAx) + T R B Tr ( x ) ) (5) 

and compare our results for each function with those of GS as well as RSE. While CE and GS [21 
differ significantly over Bc A and even more so over Bc F1 they agree with each other and with RSE [10 
on the value of Bt„- It was also only for Bc A and Bc F that RSE reported difficulty in the fit to 
leading logarithms |17j] . This strongly suggests that the source of the disagreement lies outside of the 
calculation of Bx R ^ad is most severely manifest in that of Bc F - 

We believe that the source of disagreement is the regulation of double (i.e., soft and collinear) 
infrared singularities. Calculation of Bt r involves no such regulation since the four-fermion states 
have no soft singularities, while unique to Bc F are "ladder diagram" contributions that produce the 
double singularities least controlled by energy weighting. 

To deal with infrared singularities, the exact perturbative integrands are simplified in such a way 
as to be analytically integrable in the presence of an infrared regulator (e.g., 4 — 2e dimensions) while 
producing integrated expressions that display the same singular dependence on the regulator (e.g., 
poles in e) as do integrals of the exact integrands. The simplified integrands are also used in numerical 
integrations where the regulator is necessarily removed (e — > 0) before integration. Any such algorithm 
guarantees that the singular parts of the dependence on the regulator will be correctly calculated. 

We have found that simplifications of integrands involving double poles can produce non-singular 
(0(e°)) errors from inexact treatment of 0(1/ e) shoulders of the 0(1/ e 2 ) double poles multiplying 
terms of 0(e). Since energy weighting can reposition these shoulders in a complicated way, simplified 
EEC integrands may be especially prone to such errors. These errors cannot be corrected in any 
numerical integrals where e — > prior to integration. The subtraction method prescribes addition and 
subtraction of the same quantity but the added quantities are integrated analytically while subtracted 
quantities must be integrated numerically to cancel poles in the exact four-parton integrands. Thus 
the added and subtracted quantities may differ due to necessarily different regulation methods for 
the numerical and analytic integrals. In such cases, integration of the difference between simplified 
and exact integrands is not uniformly convergent near double poles and the integrals are finite only 
in the sense of a numerically computed average. This average will generally not be the correct result 
obtained by analytically setting e — ► after completing integration rather than before. 

As a test for these errors in our calculation, the cancellation of double singularities was examined. 
Since analytic work is difficult for the four-parton states, we have focused on tests of numerical 
convergence. The scale of the independent variable controlling the singularities was magnified by a 
factor of 10 4 in a search for instabilities and neighborhoods of double poles were divided into separately 
integrated patches to isolate divergences. While further study is required, neither test produced signs 
of non-uniform convergence or error. 

Ultimately theory must be compared with experiment, and fits of our calculation to data from 
SLD have been performed p2| . Using the procedure adopted in 0, values for a s (Mz) were derived 
using the EEC as well as the asymmetry of the EEC or AEEC: 

AEEC(x) = EEC(tt - x) ~ EEC( X )- 

Renormalization scales used were in the range 

0.0035 (EEC) \ < (_|M < 4j 

- 'total J 



0.09 (AEEC) J " \E 



and while fits using KN and CE were found to have similar /i dependence, EEC fits using the larger 
CE values for B yield a s (Mz) values smaller by about 0.005 [p2l p3[. Although all B calculations 
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yield larger a s (Mz) values from EEC fits than from AEEC fit s W , it is interesting to note that the 
two differ by 0.012 for KN, as opposed to only 0.006 for CE g, g| ||: 



ct S /cE) \Mz) — 0.118 ± Q.013(scaie) ± Q.002(hadronization) ± Q.0Q3(experiment), 
a s[cE) C) ( M z) = 0.112 ± Q.003(scale) ±0.002(hadronization) ±0.003(experiment). (6) 

While the improved agreement does not constitute evidence that our calculation is correct, it is an 
attractive and suggestive feature of the results. 

We conclude that the disagreement over the next-to-leading order contribution to the EEC has 
not been resolved. Comparison of our calculation with all that is known about the EEC shows it 
to be reasonable and numerically reliable despite disagreement with previous calculations. A more 
intensive investigation of the cancellation of double singularities combined with a possible extension 
of our knowledge of the leading logarithm expansion is needed to resolve the differences. 
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